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■ There are both theoretical and experimental uncertainties in using data 



from cosmic-ray air showers to estimate hadronic cross sections above accel- 
erator energies. We outline these problems and compare the physics used to 



temporary models. We conclude that the published cosmic-ray cross section 



^ ' extract 0"*°* from air shower data to the widely used parameterization of the 

proton-proton cross section of Donnachie and Landshoff p] and other con- 

O 

oo 

Q> ■ values do not strongly constrain 0"*° fits from lower energy accelerator data. 

a: 

R*'- INTRODUCTION 



New and proposed experiments to study the cosmic-ray spectrum up to 10 20 eV and 
beyond 0-0] will depend for their interpretation on extrapolations of models of hadronic 
interactions more than two orders of magnitude in center of mass energy beyond what is 
accessible with present colliders. The interaction lengths of hadrons in the atmosphere, and 
hence their cross sections, are the most obvious determining factor in the rate at which the 
showers develop. An extra source of model dependence is the relation between hadron cross 
sections in air and the more basic hadron-hadron cross sections. 
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Cosmic-ray measurements have been used in the past to determine <r^ air and, with the 
help of Glauber multiple scattering theory 0, to estimate a^. Frequently quoted examples 
are the Fly's Eye experiment PJT0|] and the Akeno experiment ||11|| . Both experiments find 
rather large central values of <J^ aiT (~ 540 mb || and « 570 mb [|11|] at lab energy Eq ~ 
4 x10 s GeV). In both experiments, the proton-air cross section has to be inferred from some 
measure of the attenuation of the rate of showers deep in the atmosphere. The measured 
attenuation depends on the cross section which determines the depth at which showers are 
initiated, but it also depends very significantly on the rate at which energy is dissipated in 
the subsequent atmospheric cascades. For this reason, a simulation which includes a full 
representation of the hadronic interactions in the cascade is needed. Because these two 
experiments measure the attenuation in quite different ways, the fact that their inferred 
values of <r^ air agree is a non-trivial result. 

Having determined er^ ir , the experimental groups go on to derive corresponding values 
for oJJ of 120 mb |@jnj and 125 mb (DJ at y/s about 30 TeV. As noted in the Review of 



Particle Physics ||12|| , o~*° ~ 120 mb is in good agreement with extrapolation of the parame- 



terization of Donnachie and Landshoff (DL) jjj]. As we discuss in the next section, however, 
the cosmic-ray values of cr^f are based on a parameterization of the nucleon-nucleon scatter- 
ing amplitude that is in disagreement with high energy collider data. Therefore, the quoted 
values cannot be used to pin down a high energy extrapolation of the pp cross section. 

Indeed, it has been pointed out in the past [0,0] that such large values of cr™ l air 
(~ 550 mb) would require significantly larger values of 0"*°* than that predicted by the 
parameterization of Ref. [I]. Conversely, if that predicted behavior of the hadronic cross 
section is correct, then the hadron-air cross sections should be smaller, and this could have 
important consequences for development of high energy cascades. 

The plan of the paper is as follows. In the next section we discuss the relation between 
the nucleon-nucleon cross section and the nucleon-nucleus cross section, in particular, how it 
depends on the slope of the elastic pp cross section. Next we review how the hadron-air cross 
sections are inferred from air shower experiments and discuss the resulting uncertainties in 



al p-air an d their implications for 0"*°*. 



pp 



PROTON-PROTON VS. PROTON AIR CROSS SECTION 



The relation between the hadron-nucleon cross section and the corresponding hadron- 
nucleus cross section depends significantly on the elastic slope parameter B(s) 



B(s) 



d_ 

dt 



dt 



(1) 



J t=o 



This relation is discussed in the context of cosmic-ray cascades in detail in Ref . |~3| . Qualita- 



tively, the relation is such that for a given value of a ° , a larger value of the slope parameter 



inel 
air ' 



corresponds to a larger proton-air cross section. Conversely for a given value of a 1 * 
a larger value of B(s) leads to a smaller value of o*££. In addition, the smaller the slope 
parameter, the larger is the uncertainty in the derived proton-proton cross section. 

For example, the Fly's Eye value of <r™ = 122 ± 11 mb at ^ = 30 TeV 0,0 is obtained 



using a outdated geometrical scaling fit |~5],[T6| to extrapolate the slope parameter to this 
energy. This results in a large value of B > 30 GeV -2 and hence (for a measured value of 
a p n ~air ~ 540 ± 50 mb) a small value of er*"*. Using a different model for the slope parameter 
17 , 18] , for example, as advocated in the review article of Block and Cahn ||19|| , leads to a 



slower increase in B(s) and to a considerably larger value of a, 



175l|g m b El- The 



same applies to the Akeno analysis and numbers 11 



Before discussing the slope parameter further, it is useful to review briefly the basis of 
the very successful DL fits of cross sections, which are based on a one-pomeron exchange 



model (e.g. [0 and Refs. therein). In such a model, the energy dependence of the total 
cross section for AB scattering is given by |jl| 

<B(s)=X AB (^) A + Y A J^y e . (2) 
\S J \s / 

The constants X A b and Yab are target and projectile specific whereas the effective powers 
A m 0.08 and e ~ 0.45 are independent of the considered particles A and B. Within 
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the uncertainties of the measurements, this parameterization is in agreement with almost 
all currently available data on pp, pp, up, ^p, and 77 total cross sections. It should be 
noted that the high energy pp data are not fully self-consistent. There is some disagreement 
between measurements of the total cross section at y/s = 1800 GeV. Whereas the E710 



2l| and the preliminary E811 |22|] data are in perfect agreement with the DL prediction 



IJ, the CDF measurement |23[ shows a steeper rise of the total pp cross section. New data 
from HERA (a^, [0]) and LEP2 (cr*^, pjfl), although being compatible with an energy 
dependence of A a; 0.08, indicate that the cross section may rise faster with energy than 



assumed in the DL fit. Furthermore, in a recent fit to pp and pp data |26[ a slightly higher 
value of A = 0.096±g;g^| was found. 

Given the success of the one-pomeron exchange model in predicting the total cross sec- 
tion, one might apply it to derive further predictions. The one-pomeron amplitude can be 
written as 

A(s,t)=g AB {t)(-) (3) 

XSq/ 

with a(t = 0) = 1 + A. Collider data on elastic scattering suggest for small \t\ the functional 
dependence gABif) = Xab exp{|-Bo^}- Following the predictions of Regge theory, B is an 
energy-independent constant. Consequently, the elastic slope B(s) is given by 

B(s) = B + 2a'{0)ln(-) (4) 

\Sq/ 

where the parameter a'(0) is a constant and has to be determined from data [pO|j . The 
elastic cross section follows from 

( tot \2 

(1 , „2\ \ U AB) (r\ 

aAB = {1 + p) i^W)- (5) 

At high energies the ratio p between the real and the imaginary part of the forward scattering 
amplitude is small and p 2 can be neglected. 

In a model with geometrical scaling it is assumed that the increase of the total cross 
section stems entirely from an increase of the transverse size of the scattering particles. The 



opacity of the particles is considered as constant. A direct consequence of this assumption 
is the energy independence of the ratio R = crfJs) / 'a^ '(s) , which, in combination with 



Eq. (pi), leads to the relation 



167T R 

Over the ISR energy range R ss 0.17, which was the value used in @ in Refs. [P,[lT|. 
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FIG. 1. Data on pp and pp interactions JT|/ are compared with the DL parameterization Jl]/ 
(lower curve) and the fit of Ref. %2(\] (upper curve). The predictions for the elastic cross section 
from Eq. (fjj and in the case of geometrical scaling (dotted curve) are also shown. The data point 
at yfs = 30 TeV is the original Fly 's Eye estimate j^/. 



In Fig. |l| the parameterizations of Refs. Jl| and [^6[ are compared to data. The data point 
at a/s = 30 TeV is the original Fly's Eye estimate 0. The prediction for geometrical scaling 
has been calculated using the DL model for the total cross section. Whereas both Regge 
parameterizations are in agreement with data on total as well as elastic cross sections, the 
geometrical scaling model fails to describe the elastic scattering data. This becomes even 
more obvious if one considers the predictions for the energy dependence of the elastic slope 



parameter as shown in Fig. ^|. In contrast, the single pomeron exchange model is in very good 
agreement with collider data. Such an a+b ln(s) extrapolation of the slope parameter is often 
used to fit data (for example, |27[|) and also to estimate cross sections and interaction lengths 
for cascade calculations PB[p5|| . Remarkably, the minijet calculation of Block, Halzen and 
Margolis |St| (BHM) predicts a slope parameter that almost coincides with the one-pomeron 
model extrapolation using a'(0) = 0.3 GeV~ 2 . 
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FIG. 2. Elastic slope parameter for pp and pp interactions. The solid lines are the predictions 
of the one-pomeron exchange model with a'(0) = 0.25 and 0.3 GeV~ 2 . The dotted line corresponds 
to geometrical scaling. The data are taken from Refs. IS^-lMj. 



As recognized by DL the single pomeron exchange model is not consistent with unitarity. 
One way to see this is to note from Eqs. (Q,f||)D, that at asymptotically high energy the 
unitarity requirement 



a 



tot ^ o 
AB z 



(7) 



is violated. We point out, however, that the model of BHM [gOJ does satisfy unitarity and 
it gives a similar prediction to the single pomeron fit over the energy range shown in Fig. 
We summarize some of the results of this section in Fig. [3], by displaying them in the 
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( a pp~B) plane. The shaded region corresponds to the region excluded by the unitarity 
constraint of Eq. ([?[). The points represent experimental measurements at ISR (triangles) 
and pp collider (squares). The dotted line indicates the relation between B and (7*°* predicted 
by geometrical scaling with R = 0.17 in Eq. (0). This line fails to describe the highest energy 
measurements. The dashed line corresponds to the DL fit to together with equation (|J) 
for the energy dependence of the slope (with a'(0) = 0.3 GeV~ 2 ). Each point on the dashed 
line corresponds to a value of the center of mass energy of the pp (or pp) reaction. We have 
indicated with a circle the point for \/s = 30 TeV. 
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FIG. 3. B dependence on 0"*°* and the values of allowed by the Fly's Eye measurement. The 
shaded area is excluded by the unitarity constraint. Solid symbols give experimental data points. 
Dashed line shows B as in DL fit; dotted line shows geometrical scaling. The open point indicates 
°~pp °^ V$ = 30 TeV from the DL fit. The five curved lines show the region allowed by o^°^ ir = 
540 mb ±i<r and ±2a (see text). 

Using the Glauber formalism a fixed value of the p-air cross section can be represented 
as a curve in the (a^-B) plane. The five curved lines in Fig. |3] indicate the set of values of 
(Tpp* and B that result in a proton-air cross section of cr^fj^ of 540, 540 ± 50 and 540 ± 100 
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mb, that is the central value and ±1,2 standard deviations of the Fly's Eye measurement 
at y/s = 30 TeV. The intersections of the curves corresponding to 590 and 490 mb with 
the dotted line that describes geometrical scaling give the (one standard deviation) allowed 
interval for 0"*°*, as estimated in the original Fly's Eye publication. However it is clear that 
any reasonable extrapolation of the collider data (for the B—a%£) will result in the estimate 
of a higher central value for the pp cross and in larger uncertainty. Nominally the prediction 
of Donnachie and Landshoff for 0"*°* at \fs = 30 TeV is one standard deviation below the 
Fly's Eye measurement. 

It is important to notice that the experimentally measured and published inelastic p-air 
cross section is only that part of the total cross section which belongs to particle production. 
Following [13] we write this cross section as 



prod _ tot _ el _ .-.Q— e ' (Q\ 
® p—air ® p—air ® p—air ^ p—air j v"/ 

where <JpZ e a \ r is the quasielastic p-air cross section corresponding to scattering processes 
where the nucleus gets excited without direct particle production. The Glauber formalism || 
gives explicit expressions for all terms in Eq. ||. Unfortunately, there is ambiguity in the 
literature about the designation of the production cross section. It has also been called 
a p n -\ir m experimental P-[TT|] and theoretical |13| papers and it is also often referred to as 



absorptive cross section |3J-[36],[T4| . In the hope of removing this confusion, we introduce 
the notation 'prod' to represent the inelastic cross section in which at least one new hadron 
is produced in addition to nuclear fragments. 

UNCERTAINTIES IN THE P AIR CROSS SECTION MEASUREMENT 

In addition to uncertainties in converting from (7^°^ to (7*°*, there are significant uncer- 



tainties in the determination of cr^°^ ir itself. Both at Fly's Eye || and at Akeno ||11|| , the 
approach is to look at the frequency of deeply penetrating showers and to assign a corre- 
sponding attenuation length (A) on the assumption that, for a given energy, the most deeply 
penetrating showers are initiated by protons. 

8 



The Fly's Eye group measures the depth of maximum development (AT max ) distribution 
for air showers in a relatively narrow interval of S mauX , where S max oc E is the shower size 
at maximum. The tail of that distribution, well after its peak, is a measure of the depth of 
the first interaction convoluted with the intrinsic fluctuations in the shower development. 

The Akeno group selects deeply penetrating showers by cutting on showers with the 
highest size, S, at the observation level in narrow bins of the shower muon size S^. The 
reason for this procedure is that is nearly proportional to the primary energy Eq. A is then 
derived from the frequency of such showers at different zenith angles, i.e. from the decrease 
of the frequency with atmospheric depth, which is a different measure of attenuation from 
that used in the Fly's Eye approach. 

The model-dependence then is compressed into a single parameter a > 1 in the relation 

14 5 Tfi 

A = a x \ p _ air = a x prod p . (9) 

® p—air 

Here \- a ir is the interaction length of protons in air, which has a mean atomic mass of 
14.5. The effective value of a for proton initiated showers depends on the pion inelastic 
cross section in air and on the inclusive cross sections in the proton and pion inelastic 
interactions P7| , |5E|] . 

The Fly's Eye proton air cross section value of 540±50 mb is derived by fitting the tail 
of the X max distribution to an exponential with a slope of A = 70±6 g/cm 2 and then using 

are 



1.60, which is similar to the value calculated in Ref. 39 . The A values in Ref. [39 



calculated by simulating air showers assuming different energy dependences of cr^air an d 



fitting the tails of the resulting X max distributions. A values are then compared to d^°air 
at primary energy E >3xl0 17 eV. The calculation was performed with an essentially pp 
scaling interaction model ||40|| . 

Models with even very modest scaling violation, that also account for the nuclear target 
effect yield smaller values of a. It is not possible to separate the effects of the energy depen- 
dence of the inelastic cross section from those of the scaling violation in the 'one parameter' 
approach. The relevant parameter in p-air interactions is the rate of energy dissipation by 
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the primary proton K™ l air / X p _ air . The inelasticity coefficient K p n 3 air = E °~^ L ^ , where E 
is the primary proton energy in the lab system and {El) is the average lab system energy of 
the leading nucleon. The equally strong contribution of 7r-air collisions is even more difficult 
to quantify in simple terms. On the other hand a tends to saturate for very strong scaling 
violation models, because the nuclear target effects in such models are small. 



The Akeno experiment uses calculations [[11]] made also with a model implementing radial 
scaling. In Ref. |38| the results of Akeno have been reanalyzed making use of an interaction 
model with scaling violations, resulting in the derivation of lower values for a that used by 
the Akeno experiment. 

We have updated the calculations of Ref. JlJj to illustrate how different values of of^ r 
can be extracted from the same measured value of A depending on the inclusive cross sections 
of the interaction model. The calculations were performed with three interaction models 
characterized by Kf3 air : the scaling model of Hillas g@, SIBYLL g| and a SIBYLL-based 
model with significantly stronger scaling violation in pp interactions (High-K). All three 
calculations use the same input cr¥f°air = 520 mb (X p - a ir = 46g/cm 2 ) at y/s = 30 TeV. The 
resulting values of a are given in Table |. The last column of the table gives the values of 
(Jp^air that would be inferred from the Fly's Eye measurements if the corresponding value 
of a had been used. The effects of scaling violation on the shower attenuation rate used 



by the Akeno experiment [JTTJ are similar, although the numerical values of a are somewhat 
different. 

TABLE I. Cross section values that can be extracted from the measured A = 70±6 g/cm 2 with 
different interaction models. 



Model 



{K p f) 



{^p—air I 



a(Vs = 30 TeV) 



prod l 
CT p-air> mb 



Hillas 

SIBYLL 

High-K 



0.57 
0.64 



0.50 
0.67 
0.74 



1.47±0.05 
1.20±0.05 
1.12±0.05 



504 
411 

384 
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DISCUSSION 



Cosmic-ray experiments detect air showers that result from interactions of particles 
with energy up to and exceeding 10 11 GeV. Such observations have the potential to provide 
information about the growth of (T^° air up to y/s ~ 10 5 GeV. The long lever arm would be 
helpful for discriminating among models that give nearly identical results at lower energy. 
Here we attempt to summarize the problems and complications involved in the measurement 
and interpretation of cr^°^ ir in cosmic ray experiments. 

The experimental shower sets are inevitably contaminated by showers initiated by heavier 
nuclei. Neglecting this contamination would result in an overestimate of cr£i°^ r . To minimize 
this contamination, the Fly's Eye cross section was estimated by analyzing only the most 
penetrating showers, that is a subset of 20% of the entire data sample, strongly enriched in 
protons. A subsequent analysis |43|] found that the composition of primary cosmic rays may 
be very heavy in the energy region considered. If so, the contamination of heavy primaries 
could be larger than what was estimated in the original work leading to an overestimate of 

prod 
^p—air' 

The cross section estimates in Ref. |9HTT| were based on interaction models with scaling 
particle momentum distributions. Models with scaling violations predict faster shower de- 
velopment (e.g. smaller values of a). If such models were used they would imply a smaller 
p-air cross section (as illustrated in Table |). In addition, such models could also be con- 
sistent with a smaller fraction of heavy nuclei. If the shower development is described with 
a single parameter, as done in the first generation cross section estimates, it is impossible 
to distinguish between the effects of the proton and pion cross sections and the inclusive 
distributions of the secondary particles. 

Once <7p T °air is determined, the Glauber formalism can be used to infer cr^f with ex- 
trapolations for B(s) based on all available collider data. Previous analyses P-plI used a 
parameterization based on data up through ISR energies which fails to describe recent high 



energy measurements and leads to an underestimation of a, 



tot 

pp ' 
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Our basic conclusion is that cosmic-ray values of cx^f do not at present strongly constrain 
extrapolations of fits of this cross section up to collider energies. With the prospect of much 
more precise experimental measurements forthcoming from the high-resolution Fly's Eye 



and other proposed experiments [P| there is the potential for much better estimates of the 
proton-proton cross section. Realizing this potential will depend also on the use of a new 
generation [|45| of shower simulations based on interaction models that incorporate all the 
physics of minimum bias interactions up to collider energies and a correspondingly detailed 
treatment of nuclear effects. The corresponding analysis should involve a full Monte Carlo 
simulation of each experimental data set rather than characterizing the simulation with a 
single parameter. 
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